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I. INTRODUCTION 

The phenomenology of light flavor pseudcocalar mescos rj and rf provides a valuable window on many important 
non perturbative features of Quantum Chroznodynamics (QCD). It includes such important aspects as: 

• The spontaneous breaking of chiral symmetry, which gives rise to the appearance of the multiple* of light 
pseudoscalar mesons. 

• The U (1)* anomaly of strong inter actions, which gives m a s s to the singlet rfj in Nc = 3 QCD, even in the chiral 


• The explicit S(f(3y flavor symmetry breaking, due to the splitting m t / m between the strange and up/down 
quark masses (the isospin limit, where m u = m* — m and the electromagnetic corrections are neglected, will be 
assumed all through the article). 

• The 1 /N C expansion of QCD in the limit of large iVc. with Nc the number of colors in QCD. 

The interaction between the pseudoNambu-Goldstonc bcsocs (pNGBs) (r, A. rp») from the spontaneous chiral 
symmetry breaking can be systematically described through a low-energy effective field theory (EFT) based on 
5C/(3)i x 5C/(3)jr chiral symmetry namely Chiral Perturbation Theory (*PT) (1). Following large-Ac arguments [2], 
this approach was later extended, incorporating the singlet into a 1/(3) yPT Lagrangian [3—S). This combination 
erf *PT and the If No expansion provides a consistent framework which addresses all the previous issues. 

More precisely, in this article we show that this large- N c \PT framework yields an excellent description of the q and 
rf masses from lattice simulations at different light-quark masses [9-13). Constraints from phenomenological studies 
of p, a/. £ J/t> decays [14, IS) and kaon mass lattice s imu lations [16. 17) arc compatible and easily accommodated 
in a joint fit. The problems arise when one trim to also describe lattice simulations for F w% Fk and Fk/F t [10-18). 
Nevertheless, the issue of these observables in yPT is known and has been widely discussed in previous bibliogra¬ 
phy [19-23). It constitutes a problem in its own and it b not the central goal of this article. It is discussed fee sake 
of completeness and to show its impact in a global fit. 

The t) and rf mesons not only attract much attention from the chiral community but they have been also intensively 
scrutinized in lattice QCD sinmlations. where enormous progresses have been recently made by different groups [9- 
13). Varying the light-quark mooses m and m - . both their ma s s es and mixing angles have been extracted in the 
range 2C0 MeV < m, <700 MeV. We will focus on the simulation points with m r < 500 MeV in the present 
week. By observing the dependence of these observables with the light~quark mass es we will determine the \PT low 
energy constants (LECs) and further constrain the theoretical models. At the practical level we haw recast all m 



dependencies in terms of m r and study the ofcscrvabJra os functions of m*. The q and rf lattice simulations have 
not been thoroughly analyzed in the chiral framework yet and it b the central goal of the present work. Honevv, 
the numerical uncertainties resulting from our analysers in this work must be taken with a grain of salt as correlations 
between the different lattice data points and other systematic errors are not considered here. 

In addition to lattice QCD. there are also pbenccnenologxal studies of the q and rf mixing, which has been 
investigated in radiative decays <c£ light-flavor vector reso n a n ces p,u>. <S and J/tV -* VP % Py \wocesses [14. 
In these works, the modern tw\> mixing-angle scheme for the q and rf msons. which was first advocated 
in Refs. [4. 5|. was employed to fit various experimental data. The common methodology in these works is that 
the two-mixing-angle pattern for the q and rf b simply adopted to perform the phenomenological discussion and 
the mixing parameters arc then directly determined from data. Thb b a bottom-up approach to address the q-rf 
mixing problem and it b quite useful for the phenomenological analysis. Contrary to the bottom-up method, it b also 
very interesting to study the 17 - 17 ’ mixing from a top-down approach in which one first constructs the relevant yPT 
Lagrangian and then calculates the q-rf mixing pattern and parameters in terms of the LECs. In thb case, one con 
predict the ryq* mixing parameters coce the values of the unknown LECs are given. The present work belongs to the 
latter category of top-down approaches. 

Though the singlet meson, which is the m am component of the ph>mcal q* state, b not a pNGB due to the 
strong f/(l>4 anomaly, it can be formally introduced into \PT from the large -No point of view. The argument b that 
the quark loop induced U(l) A anomaly, which b responsible for the large mans of the singlet rjj. b 1 /Sq suppressed 
and hence the f)j becomes the ninth pNGB in the large Sc limit |31|. Based on thb argument, the leading-order 
(LO) effective Lagrangian for t/(3) xPT. which simultaneously includes the pNGB octet ir. K. q* and the singlet rju as 
dynamical fiefds. was formulated in Ref. [3]. Later co. a full 0(p A ) U{ 3) chiral Lagrangian was constructed in Ret (7) 
and the discussion on the 0{p 6 ) unitary group chiral lagrangian has been very rececfcly completed in Ret (32). Subtle 
problems about the choice of suitable variables for the higher order U{ 3) \FT Lagrangian in the large Afc framework 
were analyzed in Ref. ( 6 ). 

The standard power counting employed in SU{ 2 ) and 51/(3) in poms of the external momenta and quark 
masses [ 1 ], b not valid any more in C/(3) \PT. due to the appearance cf the large imass. However, since the singlet 
»\j moss squared b ch iW B like if Sc in large S c limit, the rfo mass can be harmonized with the other two expansion 
parameters if one assigns the same counting to if Sc, the squared momenta p 2 and the light quark masses m,. As 
a result of thb. in order to have a systematic power counting, the combined expansions on momentum, light quark 
mass es and i/Sc arc mandatory in 1/(3) *PT | 6 . 7). We will work in this combined expansion in our study and 
denote it as 6 expansion throughout the paper, where 0(6) ~ 0{p 2 ) ~ ~ 0(1/No)- Thb counting rule b 

different from the one proposed in Ref. [33]. where the »mass b counted as 0 ( 1 ) and the infrared regularization 
method b employed to handle the chiral loops. 

Some recent works in Refs. [ 8 . 34-37) have addressed the 17 - 17 ' mixing in the chiral framework up to next-to-leading 
order (NLO). As on improvement, we will perform the systematic study of the q-rf mixing in the 6 -expansion scheme 
up to next-to-nert-to-leading order (NNLO) and take 11 x 0 account the very recent lattice simulation data, which are 
not considered in the previous weeks [ 8 , 34-37). In addition, we also simultaneously analyze the m r dependences of 
other physical quantities from lattice stnmlatioos. such as the axial r. K decay constants and the mass ratio of the 
strange and up/down quarks, in order to further constrain the yPT LECs. 

Thb article b organized as follows. In Sect. II. we introduce the theoretical framework and calculate the rcletunt 
physical quantities. In Sect. III. the phenomenological discussions will be presented. Conclusions will be given in 
Sect. IV. F\irther details about the calculations up to NNLO arc relegated to App. A. 



II. THEORETICAL FRAMEWORK 
A. Relevant chiral Lagrangian 

At leading order in the 6 expansion, i.c. 0 ( 6 °). the f/(3) xPT Lagrangian consists of three operators 

- xK.u'') * £(*-> + . <1) 

where the dural buildup block* ore defined u [ 1 . 6-8] 

x = 2B(. + .», x»=«V±«xV X ■ t>£ (det U ), 

u„ « i, D.U = 0„U - i<i'„ + a K )U + iV{v M - a„). 


(2) 
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irith the pNGB octet+singlet matrix 
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3 » p. being the external senior. peeodascalar. vector and axial-vector sources, respectively. The coupling F 

Appearing in Eqs.(l) and (2) corresponds to the pNGB axial decay constant in the large Sc and chiral limits. The 
light quark mor®« are introduced by setting (s + ip) — diag(m t m. m,}. being rh the averaged up and down quark 
ma s s es and m t that of the strange quark. 

Notice the structure of the LO Lagrangian in Eq. (1): the first operator is of 0{Sc,p l ) tjpe. the second one 
coTTfapocds to the type of 0(Nc, m f ) and the last one steins from the QCD f/(l)a anomaly and is of O(S^.p 0 ) 
type, where U is counted as 0(1),** ~ 0{Nc) And AfJ ~ 0(JV£ l ) in the classification 0(A r £,p k . mj) of the EFT 
Lagrangian operators in Eq. (1). In the following, we will denote the chiral expansions in powers of squared momenta 
p 2 and quark masses m fl simply as a generic expansion in p 2 . 

The NLO U (31 chiral Lagrangian. i.e.. 0(4). contains 0(A T c.p 4 ) and OiS^.p*) operators. The relevant ones in 
our work read [6 


c, ° - + T<x*x+ + x-x-> * ZJr-rrxD'X - ££x<x _). 


F 1 A, 


W 


with the dimrtuionlraa LECs’ * aline lilcc Li.L,~ 0(iS c ) and Aj.A* — 

At NNLO. i.e. 0(A J ). there arc three type* of operators; OfNj'.p 3 ). and OfNc-.p"). Their raqriKit 

forms read |7.18] 

+£,(u*‘U|)(x.) + I»(X4.)<X*.) + MX-){X-) + iu(“»K““X.) + ^X(x+x-) 

+Cx rlV^X.) + Ou(v*' , X + X*> + C„^ M x^X.) 

+^«(X.X.X4> + C S1<X-X-X*). <5) 

"here the first line corresponds to the 0 (Nc‘.p‘) type, the seeood line is of the ©(Aj.p 4 ) type and the last two 
lines are of the 0{Nc,p r ‘) type. The LECs carry the scalinrs Lt.Ls.Lr.Lis.Ijs — 0(N%) and 

Cir.C„.C 17 .C„.C„ — O(Nc). Notice that we haw only shown the operators at different i orders in Eqs. (I), (4) 
and (5) that arc pertinent to our present study, not aiming at giving the complete sets of operators. The conventions 
to label the IX). NLO and NNLO operators in Eqs. (1). 14) and (5) follow closely the notations in Refs. (6, 7. 38|. 
Unle» it is explicitly stated, the LECs will correspond to C/(3) *PT and must not be confused with thc«e in SU( 3) 
*PT. The matching between these two EFTs can be found in Ref. [6|. The terms L } are denoted as ^ in Ref. (7. 8|. 

Comparing the t/(3) and SU{ 3) theories one can observe that some terms have been reshuffled in the A expansion 
of the C/(3) Lagrangian. For example, the L l - 4 .s.e.Ta terms are NLO in SU{3) *PT. but they are now split into NLO 
and NNLO in the A expansion (see Exp. (4) and (5)). We have several additional new operators, namely the last one in 
Eq. (1). the A*»u in Eq. (4) and the r£ 2) .Lu,L» terms in Eq. (5), that arc absent in the SU{ 3) *PT case. Finally, 
the chiral loop* start contributing at NNLO in the A expansion, while they appear at NLO in the conventional SU( 3) 
case. 


B. The r?+( mixing at NNLO In <J expansion 


Next we calculate the v^r/ mixing order by order in the A expansion. In literature, there arc two bases to address the 
Tj-rf mixing, namely the singlct-octct basis with ffc and and the quark- flavor basis with q, and rj a . The relations 
between fields in these two bases ore 



In the large-AV limit where the U(1)a anomaly is absent. q T and rj M are the mass eigenstates and they are generated 
by the axial-vector currents with the quark flavors qq m (uti + dd)/v^I and sf, respectively. The two bases are related 
to each ocher through an orthogonal transformation and provide an equhnlent description for the rf-rf mixing. 



As noticed in Refs. ( 37 . 39 j. when doing; the loop calculations with 7 and rf , it is rather cumbersome to work with 
the 73 and 7* states. The reason is that at leading order the Lagrangian in Eq. (1) gives the mixing between 73 and 
T} 4% and the mixing strength is proportional to - m 2 . which in the 6 expansion is formally counted as the some 
order as the diagonal terms in the mass matrix for 73 and 7*. As a result, the insertion of the 73-7* mixing in the 
chiral loops will nee increase the i 5 order of the loop diagrams. This makes the loop calculation technically much 
more complicated, as one needs to consider the arbitrary insertions of the 73*71 mixing in the chiral loop diagrams. 
Nevertheless, Refs. | 37 . 39 ] provide a simple recipe to handle this problem by expresssing the Lagrangian in terms of 
the Tf and T f states which result from the diagooalixatxm of rjo and 7* at leading order in < 5 . The main difference is 
that the mixing between 1 J and if is now at least a NLO effect in 4 . while the 73*7* mixing was appearing at LO. The 
relation between the IX) mass eigenstates 7 and rf and the singlet-octet basis is given by the mixin g angle 0: 



( 7 ) 


with c# = cos $ and s* = sin #. The LO mixing angle 0 and masses of 1 J and if are given by the leading order 
U?nmB*n in Eq. (1) (*« e.g. Ref. | 37 j): 
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with A 2 — Tlfjf - TIT 2 . Here TTT/c and TIT, denote the LO kaoa and peon mooses, respectively. 

When higher order corrections are taken into account, the IX) diagcoolized 1 J and if will get mixed 
the NNLO, a general porametrixation of the bilinear terms involving the 17 and Tf states can be written 

c = i ±d M if^~^'n*Y d ~ d ‘ y,d ~ iryr * s> 


«> 


l-t-V 




1 + *I' 




mi-f J 


Sn,- 


"4 + 4 -t 


TV - W - 


(ii) 


where the tfa contain the NLO and NNLO corrections. Here these operators must be understood as the terms of 
the effective action that provide the pseudeoealar meson self-energies. The higher-derivative terms the 

first line of Eq. (11) are exclusively contributed by the 0 {p*) operatcc C i2 in Eq. ( 5 ). which belongs to the NNLO 
Lagrangian. The remaining %a receive contributions from the NLO operators in Eq. ( 4 ), the NNLO ones in Eq. ( 5 ) 
and the one-loop diagrams, which contribute at NNLO. Their explicit expressions con be found in App. A. 

At leading order, there is only the mass mixing term from Eq. (1) whereas at NLO and NNLO ooe has to deal 
in addition with the kinematic mixing terms in Eq. ( 11 ). apart from the maw mixing. The physical states of 7 and 
rf can be obtained from the perturbative-expansion (X expansion) in three steps: os a first step, we eliminate the 
higher-derivative terms through the 6eld redefinitions of rj and rf: then we transform and rescale the fields resulting 
frem the first step in order to write the kinematic terms in the canonical form: after the preceding two steps, there is 
only the mass mixing term left, which is straightforward to handle. 

In the first step, we make the following field redefinitions fee the 17 and Tf states 


T?-. 7 T+Q,Oy+Q 3 caj' < V -of *a-20ri*a,0rf. 


( 12 ) 


with the d'Alembert operator □ 2 d M d*. After some algebra manipulations. it is 



aight forward to obtain 


( 13 ) 


so that the three higher-derivative terms in Eq. (11) will be eliminated. Notice that the Oj 2 , are NNLO. i.e., 0 {P). 
Substituting the field redefinitions from Eq. (12) into the general mixing structure in Eq. (11) and keeping the terms 



up to NNLO, the resulting bilinear Logrongian reads 


C = 




1 + 6 * + mL 


-,-WV +--f (*, + T (m^.-f <)] 

" 4 + 5 -i. w . . 

-r-=W- ~ A-*W - 


(u: 


In the second step, we need to eliminate the kinematic mixing term in Eq- (14). and then to rescale the fields to 
hove them in the canonical forms. This can be done pert urbati rely. In the final step, we take care of the mas 
term. The last two steps can be achieved through the following field transformations 


I Ti\ /eastf* -sintf*\ (l + * A ** \ I V \ 

V^J-Vsin#, cm#, ^ S„ 1 + 6 C J\’f)’ 


(IS) 


with rj, rf the physical states and 

., _ «r . "fa W 4 2 *u> 

*• Y ~2~ ~s~ ~s~- 

, J* , *J, J , 1 . A^NLOit.KLO V-'U,J«NLO 

5a " •v + T l ’"*' tm * ) -5-3-• 

fc -a., 

where ^.nlo. «V,nlo. ^ki<LO stand for the NLQ parts of the three quantities respectively. We point out that 6^ 6* 
receive both NLO and NNLO contributions, while ^,4,,6, ore only contributed by the NNLO effect, which is the C l2 
operator in Eq. (5). Comparing with the NLO results in Eq. (15) from our previous paper (37]. we hove generalized 
the expression to the NNLO cose in the promt Eq. (15). Another way to treat the mixing of pseudoscalar mesons in 
*PT was also previously studied in Ref. (40) and applied to the x°*i) case up to the two^loop 1ml. 

In the practical calculation, it is more often to use the inverse of the relations in Eq. (15). where the perturbative 
expansion leads to 



6 f b \ / castf* sin** \( rt\ 

1 +*cj V-™ 19 * CM#, ) \r/ ) ’ 


(17) 


•nth 


_ At _ "fa 1 . **n»LO . tt2.nu. 


_ «** it, 1 . 2 . . 3Ar.KLoJ.SU) 3V.NLOJ..KLO 

* B = -y'T 1 ^^) -i- + -3- 
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The #, appearing in E<p. (15) and (17) is dMcrmincd through 

t*n#i 


m 




(18) 

(19) 
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Frith 


- j4».NLO (iiniNLO •* i ^,.NLo) + g4i.NLO&r.NL0 ( 3m J-* S'"*) 

- ULO (Ar.NLO * "V.NLO> ■ 

mj = '4* 3 -,i-n4(* r -*>T^ 1 ) + n44 tNL0 -*^rT^? >Il0 -t-^mi.A; NU( 

_3 *.NL0 3 m‘.NLO “ *T.NUl 3 miNIO - 

m\ = - mji (V * tr$ Ji) + ^^.nio + {"#, 3 nlo + 

-A,.NLO<*m>.NLO - A*>U>A„*,.NLO - 


K “ *1**1 ~j (mj - <) + 4ft,. 

2m? 


"" 


"»* + if ("i " ™«) 1 * *%>• - 


( 20 ) 


where 6 ,.nlo stand for the NLO parts of V 

In the phenomenological discissions. the popular two-mixing-angle parametnzaticci in the singlet octet basis [i. 5) 

the form 


t) \ 1 ( Fn cos0* -Fosinfl* 

i I T \ faring Focos^o 


( 21 ) 


Combining Eqs. (7) and (15). it is straightforward to derive the relations between the four parameters in the 
two mixing-angle scheme m Eq. ( 21 ) and the vPT LECs: 



ton#, - 

tan 0 D = 


F 2 1 [cos(fl + 0*) + 6iisin<0 -0j) + 6a costfeosty - 5cnn^5inPj] J 
+ [n»14 + 9 j) + Aflao(» - 9 4 ) + iAcatfitwfi, + Afsm#aoS 4 ) 3 

t 

F 2 | [ - sinffl + 0s) + 6s oas(& - 0s) - d^stn^cos^ - 6ccos0sin0*) 2 
+ [«*!* + 9.) - A o am(0 - 8 t ) - S A iOx8f.ax9 t + JfCoiflao^] 3 

t 

sin(fl + flj) + 6s cos(ft - 0s) + 6 A cos05zn0s + fcsinflcqsfr 
coofJ + 0 |) + 6s *n(0 - 0s) + (xs 0 cos- 6csixi0*n5J ’ 

- sjn(fl + 0 j) + JjfCos<^ - 0 s) - 6,4 sin# cos 0 j - 6c cos0 tan 0s 
co3(6-0s) - Sbmd{0 - tfj) - 6,4 sin 0 sin Ot + 6 <.cos 0 cos 0 j * 


( 22 ) 


where the *PT LECs arc implicitly included m 0 . 0 *. 6a. 6b and 6c. Since 0 s.6a.6b.6c ~ 0(6) or 0 ( 6 2 ). at LO one 
has F 6 * F 0 = F and one mixing-angle 0* - 0 O - 0. 

The relations between the physical q. if' states and the quark-flavor basis is commonly parametrized as 



F, cm 


-F, sintf 



(2J) 



Comb min* Hoi. (6). (7) nod (IS), it e strnirfrtlorwnrd to obtain the potametera in Eq. (23): 


^ - 
f; = 

Inn 4 , = 
tnnt*, = 


2 F* * FI - 2ZlFeF, tnn(fl 0 - #,) 

-3- 

« + 2J’? + 2v*F 0 F..ta(# t -O 


v ^F,<cotl,-f fygnt> 0 
y'lFo »in Up - F, cos 0, ' 
y/lFp ccc fl D + F» nn fl» 
v/lF* sin 0 | - Fo cos ' 


(24) 

where at LO in the 6-expansion ooe has F, = F. = F and m 0^-0. with the ideal mixing Aid = - orcsin \/2/3. 


C. Insights into previous studies of the tyr/ mixing 


In the previous subsection we have performed the full computation of the mixing up to NNLO in the 6 expansion. It 
is interesting to make a brief s umm a r y of the assumptions made in previous works, where plenty of mixing formalisms 
have been proposed to address the r^rf system [8. 24, 34-37, 41. 42). In Ref. |41), only the lowest order in the quark 
mosses and 1 / Sc % i*. the IX) contributions in the 6 expansion, were token into account. Even though it provided 
a reasonable first approximation, it failed to give an accurate description of the experimentally observed mo*a ratio 
mj/m 2 . The Oip 2 ) contributions were studied up to NLO in l/Nc in Ref. |42j (including the terms in Eq. (1) and Aj 
and Aj in Eq. (4)), perfectly explaining the experimental value of However, it turned out to be inadequate 

to give a proper value for the q-q* mixing angle. On the other hand, the authors in Refs. (34-36) went up to NLO in 
the p 2 expansion but keeping just the LO in 1 /Sc (including the terms in Eq. (1) and and U in Eq. (4)). Both the 
q-q* mixing angle and the ratio Fk/F V were qualitatively reproduced in this case. The full set of NLO contributions 
in the ^-expansion (i.e.. the effects up to NLO both in 1/JVc and p 3 ) was analyzed in Ref. (8). together with the 
mixing angle and the x % K. q and q' axial-vector decay constants. In Ret [37], the contributions from the tree-level 
reso n a nc e exchanges and partial NNLO effects, e.g. the loop diagrams, were considered for the masses of q and q\ In 
this work, we generalize the discussions up to the full NNLO study in the ^-expansion and confront our theoretical 
expressions with the very recent lattice simulation data and the phenomenological inputs from the two-mixing-angle 
scheme. 

Reference |24| introduced a quark-model inspired approach to the q-q' mixing, which is commonly referred as the 
FKS formalism and used in many phenomenological analyses [43]. The essence of the FKS formalism is the assumption 
that the axia l decay constants in the quark-flavor basis takes the same mixing pattern as the states 



COSO -sin* 
sin* cos * 





(25) 


where the decay constants are defined as the matrix elements erf the axial currents 

(0|^(0>|P(t» - 1 V 2 FIK . (a = 9.P.,,„'), 

K ~ . A‘ m = . (26) 

From another point of view, the pattern of Eq. (25) employed in the FKS formalism relies on the assumption that 
there is no mixing between the decay constants of the flavor states and q,. In the *PT framework, the physical 
masses and decay constants can he obtained from the bilinear parts of nonet fields in the effective action with the 
correlation function of two axial currents. Since the correlation function is the second derivative with respect to the 
axial-vector external source and a M always appears in the Lagrangian together with the partial derivative as 
shown in Eq. (2). the absence of the mixing for the q, and q f decay constants in Eq. (25) implies that there are 
no kinematic mixing terms for the quark-flavor states q, and q, in the FKS for ma lis m In fact, the assumption in 
Ref. [34) is in accord with the FKS formalism. This can be simply demonstrated by expanding the chiral operators 
considered in Refs. [34. 35), ir. those in Eq. (1) and £&.£« in Eq. (4), up to quadratic terms in q f and q*. 1 No 


in R«f> 134 36;. r*f>«*iv«ly. Tbt A Ura m tbt pr«wut r«ftr«t>:c* 
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kinematic mixing Knns fee the and i)j fields result from these chiral operators. This also confirms the finding in 
Ref. |36| that only when tbe NLO of l/Nc operatcc is excluded the FKS formalism is recovered with their chiral 
Lagrangian calculations. 

Since general terms up to NNLO in 6 expansion arc kept in our discusiico. unlike in the previous works (8, 34- 
Y7, 41. 42] where different assumptions, such as the preference of the higher ccdcr p 2 and 1 /N c effects. are made, it 
ls important and interesting for us to justify these assumptions in later discussions. 


D. Masses and decay constants of plon and koon up to NNLO In A expansion 


The NLO expression of the pion decay constant in tbe A expansion reads 


F r = F 1+41* 



(27) 


or. up to tbe precision considered, one can also use the physical F r in tbe expression inside brackets. 

(“) 

The differences between Eqs. (27) and (28) are NNLO effects. We mention that at a given order there is always 
ambiguity in chocs ing the renormalized quantities in the higher order expressions. In contrast, there is formally no 
ambiguity in the expressions in terms of the quaefcity F. which is tbe pNGB axial decay constant in the chiral and 
large Sc limits. For ex a m ple, if we limit our analysis up to NLO. formally, it is equally good to use F r or F* in the 
denominators of tbe NLO part in Eq. (28). since the difference is beyond tbe NLO precision. A typical solution in the 
chiral study is to express the quantities, such as m r . F n , itik*Fk, in terms of the renormalized F r in the higher order 
corrections, as dcoe in the two-loop calculations in SU{ 3) *PT (44). We follow this rule throughout the current work 
to estimate the uncertainty due to the truncation of the A expansion when one works at a given order in perturbation 
theory. We mention that tbe notation of mj in the above equations stands for the renormalized pion m va squared 
and the leading order m a s s squared is denoted by TTr^. Notice tbe IX) pion m a s s squared TTf£ is tbe ooe that is linear 
in tbe quark masses. The expressions relating mj and T7T 2 will be discussed below. 

Similarly up to NNLO. we can cither use F or F r in tbe NLO and NNLO expressions fee other quantities such as 
Fk and the 4/s in Eq. (11). In tbe NNLO expressions, the difference between using F cc F, in tbe denominators is 
a next-to-next to next toleading order effect (N 3 LO). Since in this work we study lattice simulation data up to pion 
mass of 500 MeV. the convergence of tbe chiral scries is expected to be much slower than that in the physical case 
with m T - 135 MeV. Therefore it is a prion not trivial to judge whether tbe two approaches-using 1/F 2 and 1/F 2 - 
are numerically equivalent or tbe Lattice data prefer one of them. Indeed in Ret (45). it is already noticed that to use 
F or F r could cause some noticeable effects. Wc will use the difference between both approaches as an estimate of 
the truncation error at a given order in 6. 

We take the pion decay constant as an example to illustrate the differences of using F and F, in tbe higher order 
expressions. Using F in the higher order corrections, its expressico reads 



F \l * + ,24 LI - 64 UL t )< + (8C, 4 + 8C, T )^ 

Ailmj) | /4 0 ,m^) 

16r 'F J 32^F l 


The cne-point loop function Ao(m 2 ) is calculated in dimensional regularization within the M2> — 1 scheme [1| and it 
reads 



(30) 


u drew**! ib Umu 
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with the rccormalizaticQ scale p fixed at 770 MeV throughout. Using Eq. (27) to replace F by F T in tbe KLO and 
MNLO corrections, the resulting form is 


F n » F 


i + + At r r F r* * * («cm +^yjf 

mj) ^ A 0 (m^) 

T^T? WTT 


mi + 2 fn 


mi 


mi 


[31) 


In the <5 expansion, the expressions for a physical quantity with F or F r in tbe higher order chiral corrections differ 
only for tbe LtLjwi.i and L\ A ; -i.i terms, since the differences by replacing F by F r arc originated from the NLO 
expressions of F„ in Eq. (28) and we only retain terms up to NNIX) in this work. It is clear that the difference between 
Eqs. (29) and (31) is tbe L{ term. Notice that in tbe <5 expansion scheme, the terms like I*L 4 are N 3 LO and will be 
dropped throughout the article. 

The corresponding explosion for the knon decay constant when one us« F to express the NLO and NNLO 
cocxcctxicfi reads 


_ __ 4 r m* mj + 2 m 3 _ a _ _ .ml 2 m* - 2m\.m\ + m* 

fit — F 1 + iU-jO- * iL t ' F , * + (24 L\ - 64 UL,)-£ * 8C M —*- 


mi(2<-m:) . 3i4o(mJ) . 3A,(mi) . &24>(”4> . 3 4M<) 

-P- + TSfr»P' + WP ' + 12U^P + ' 


+ 8 C„ 

On tbe other band, exp re Kang the NLO and NNLO contributions in term* of F n yields 
Fk - F 


i . jr m * , jf < + 2mi 

1 + y + 4L 4 fj 






+sc,. 


m{( 2 mi - mi) 


(32) 


3A,lmj) _ 3d„(mj,) . 3.'A,(m^) 

TSPFT TGTFT v&x'-Pz T^r 7 TT~ 


(33) 


The expanded expression for the ratio of F/c/F r in terms of F up to NNLO in 6 exponioaQ. takes tbe form 



F* pi 

5A.(mj) .-Mm},) . 3<3A,(m') 

IJSrSP * dirlpi “T2&7’P" “TTiSTP”" 

When expressing the previous result in terms of F v% it reads 

1 4U mi ~ m - + 8 LI*"** ~ 2m * mi ~ &mi ♦ UUL,?Lj< 

Fn F, FT Fj 


+16 C u 


mi - m 1 




T? 


18C.T 


it^l 


- m 




5A)(mJ) ^ ^ S^Arfm*) 3»jA 0 (m>.) 

" laSiT'R 64P7T lttr'K 129!r J Fi ' 


(34) 


(») 


■hich diSW» from Eq. (34) m the L\ term. 

The p»u squared moss up to NNLO B riven 


"4- T 4 + «4* u, +*4-*“. 


(36) 
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mi m 2 Bm . (37) 

2.NLO _ 8<2£i - L,)mi 

m » “ - J72 -• V»J 

m j.NNLO = WU - LtimlVml * mj) 64<iLg - 6 L-M + 

16(2C|2 + Cm + C„ - 3C„ - 2C, l )m , > mj(cj - 2v/3r # «, + 2*})ilo(*n3) 

P * %PP 

m J(2c^ + 2 >/Icrn + »j)A)(ml) m ’,t 0 (m’) 

*-5S77T’-i 2 .’P ‘ m 

When expressing the rccormahzed m T in terms of F n . the only differences Are the L\L* and L\ terms in Eq. (39) and 
the other porta are the some as in Eq. (36) with the explicit replacement of F by J> in Eq. (27). Therefore we only 
pv« the different ports for simplicity when expressing in terms of F, and they read 


Bn 2.(i.M,L..L! = 128(4L,£.,- Lg)m« 


The mass squared for kaon up to NNLO is provided by 


mi = Pli * mi NLO * mV*’"" 


mi — B(m + m,), 

2 .KLO _ 8 ( 2 I» - Li)m% 

K ~ P • 

i.NLO - U >mi ( 2 ml + m}) 64(1? - 61,1. + 8 lJ)m 


32C„mi 32C„mi 16C ir mirn;(-2mi + m{) 16C M m*(2m* - 2mim{ ♦ m*) 

~—p -P“" *-P-"-P- 

18C 19 mJ i (2mi - 2m‘ K ml * m‘) H(3m* + mi) + 2v^,»,(-2mJ c * mj) - 4m* ^ 4 ( 11 $ 

P ' I«PP 

l-lrfml + 2 v^,« r ( 2 m' - m\) + (3m' ♦ m*)»2A a (m* ) 


When expressing Eq. (41) in terms of F,. the dUferentes »rc the L\L% and Li terms in Eq. (44) and the new 




&4tjmi(mi -> mj) ( » mj) 


Notice that the masses of pico and kaon appearing in NLO and NNLO parts in the above equations correspond to 
the renormalized quantities, instead d their LO expressions. In addition, this gives the quark mass ratio relation 
m,/m » - 1. 

When performing the chiral extrapolation of the lattice data, instead of the renormalized as in the previous 
equations, it is convenient to use the LO kaon maos squared in the higher order corrections. In this way. we do not 
need to iteratively solve Eq. (41) in order to give the value of m K for a given m T . The result in terms of TTT K in the 
NLO and NNLO expressions becomes 

1U _2 _2,La-NLO XLrt-*NLO (AM 

m K = Tir K +m K + m K , (461 



with 


2 L-* 1 -.NLO 

K 


J.U-.SNLJ 

K 


8 (2U - Ufm*, 

-PJ-' 

8(2/., - i«)mi(2mi + mi) e4(Li - 2Lt.L,)m% 32C„1»& 32C„m' > 

-P- + - F* - F* - F*~ 

16C,-mimi(-2WTi ♦ mi) 16C 14 W^(2rr} r - 2TTP K m‘ * m‘) 

*- Ti -T5- 

, 48C„mi(2m< c - 2rrimi + m*) 

F‘ 

[4(3mi^mi)- t .2v^„(-2nri. + m»)- 

192r 2 F* 

(-4^^ + 2vfer*«<2m5 f - mj) + (3m*. + mJ)jJlAo(mM 


Wh:L 


Eq. (46) in terms of F„. the dffferences ore the L*I* and Li terms in Eq. (48) and the new exp 


a.UM*'.KL,L..L2 WI R (L 5 “ 2L 8 (WT^(»JI? f - mj) 

m K = -- 


When confronting with the Lattice data, we only consider the simulated points with the ph>uicol strange-quark 
maos. ie. the lattice ensembles that when extrapolating to the physical pico masses lead simultaneously to physical 
icaon masses. In this case, we con express the LO kaco mass squared as 

- .... .... TTTi*'‘ r TTTl 


mi = + m) = BlmJ"**’ + m"‘') - Bm™ + Bm - mi , ' h ’ - 


■bar Vt‘* U ' and mj- 1 ’’ can be obtained through E«p. (36) and (41) by substituting the physical tnamri of r. K. n. r/ 
in the NLO and NNLO expressions. For T7T*. which varies in the lattice simulation, we can extract its value by using 
Eq. (36). In this case. m r in Eq. (36) takes the value from lattice simulation, and m K . m 0 . which only appear in 

the NNLO part, can be approximated by their LO expressions. 

In the above discussions, we have distinguished the situations of using 1 /F 2 and 1/F^ in the higher order corrections 
for various observables. Similarly, we can also generalise the discussions by replacing the renormalized masses (m, 
and mjc) with the LO onra (7 TT t and ITT*) in the higher order corrections. We take the observables F r and F/c as 
examples to illustrate the differences. The renormalized m r and m* have been used in Eqs. (29) and (32) fee F„ 
and Fj* with 1/F 2 in the higher order terms, respectively. After replacing m T and m* in Eqs. (29) and (32) with 
their expressions in terms of the LO mosses TTT r and TTT X through Exp. (36) and (41) respectively, the corresponding 
expressions arc found to be 


*• 1 J. 4r ^ ’’'i + 2 a nrJ 

F 1 + 4/, 


|.SC„ + 8 C„)-^. + 


K . i) Ajimi) 


F 1 16r ’F 7 UrtF* ' 


Fx = Fll » 4/, 


‘ 8 Ci- 


. 1 y , TTT * + _ 8 + 8 Ci, ~' rr ' ~ 2lTT *' m * ''"'t 

m^(2mi -m^) 3A,{mi) 3^ 0 (mi) 3^K) S^4.(m^) 

r- 12 a?? + 12 &HF* 


As m the discussion of \/F 2 versus 1/F* up to the NNLO precision, the expressions for a specific observable by using 
the renormalized mosses m T .mx and the LO TTT*. TTT* only differ in the terms like L t L }% being L t and L } the NLO 
LECs in Eq. (4). This con be clearly seen when comparing Exp. (29) and (51). Eg. the differences caused by using 
the renormalized masses and the LO coes ore the L\ and L\L% terms, apart from the explicit replacement of m* and 
m* by iir r and T7T#c respectively. Similar rules arc also applied to Eqs. (32) and (52). 

To replace m F . mx by T7T* and TFTx in the NIX) and NNLO corrections in Eq. (36). the only changes happen for 
the LiLt terms and the corresponding new expressions read 


2.\W..wr*.n.L,i,.4.;.4.; 


M U{U - 2L%)irfi 
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In principle. we should also present the results expressed with the LO masses TTT* t TTr* and the renormalized decay 
constants F w% Fk . which can be straightforwardly obtained by substituting the relations in Eqs. (36) and (41) into 
the corresponding observables. We consider the expressions given in terms of the renormalized masses and l/F 2 
as our preferred ooes in this work. Tbe reason to cbocoe the renor ma l i zed m a sse s is for practical purpose, since in 
lattice simulations tbe different observables are typically given a a functions of the renormalized mj. Also most of the 
chiral studies choose to express the quarfcitics with the renormalized masses in tbe higher order corrections, such a a 
in Refs. [44. 45]. Following this rule we consider tbe results with tbe renormalized masses and l/F? as an estimate of 
systematic errors due to the truncation of tbe 6 expansion when one works at a given order in perturbation theory. 
While for the case with tbe IX) masses, we shall also comment the results in tbe following numerical discussions. 


III. PHENOMENOLOGICAL DISCUSSIONS 

The beg challenge in tbe present general discussions on tbe q*q' mixing is the determination of tbe unknown LECs 
in Eqs. (1). (4) arid (5). The recent lattice simulations co the light pseudeoealnr mesons provide us valuable sources 
to constrain these free parameters. The considered lattice simulatxms include tbe m n dependences of tbe marges 
of |9-13j and kaon (16, 17]. and tbe x.K decay constants (16, 17] and their ratics (18]. Moreover, relevant 
phenomenological results and experimental data will be also included to constrain the LECs. 

Since we do net consider the isospin violating effects, we will take the values for the physical pion and kaon morses 
in tbe isospin limit from Ref. [46], where the corrections from the electromagnetic contributions arc removed. 

m T = 135.0 MeV f m K » 494.2 MeV. (54) 

These values will be used in later chiral extrapolations, while for the physical masses of ij and rf and the decay 
constants <X pion and knon. we will take tbeir world-average values from Ref. [47]. 

In order to show the results step by step, we present tbe discussions in the following sections split in three parts: 
we consider fits performed at leading order, next-to-leading order and next- to-next- u> leading order. 


A. Leading-order analyses 

At leading order, tbe Tf-rf mixing b described by one free parameter, namely tbe singlet qo mars Mo in Eq. (1) 
and tbe explicit expressions fee the masses and mixing angle are given in Eqs. (8), (9) and (10). At this order, the 
r. K decay constants are degenerate and given by their chiral and large Nc limits. Le. F r — F* — F. Therefore we 
shall cot take tbe lattice simulations of the decay constants into account for tbe LO discussion as they clearly show 
the need of higher order corrections for a suitable description. Also at leading order. F will not enter tbe masses and 
mixing angle, as shown in Eqs. (8). (9) and (10). As a result of this, we do cot need to distinguish tbe two situations 
with F ox F x discussed previously in tbe expressions of different observables. Apart from the lattice simulation data, 
we also fit the physical values of tbe q and q' mooses. Nonetheless, fitting the physical masses with the experimental 
precision at the level of several hundred-thousandth is too ambitious. Since tbe ultimate goal of the present work 
b the NNLO study, the ballpark estimate ci our theoretical uncertainty, starting from the N’LO part, should be 
around 3%. This value b obtained from the general rule that each higher order correction in 6 expansion, either the 
5C/(3)-flavor breaking cc tbe 1/JVfc «£<ci. b around 30%. In fact, tbe estimated three-percent uncertainty b also 
similar to the typical error bars reported in many lattice simulations, in the range from 3 - 10% (9-13]. Consistently, 
we assign a 1 % uncertainty to tbe physical values of m v and m,y in the fits. 

The value c£ the singlet mass Mo from the IX) fit b 

Afo = (£35.7 ± 75) MeV. (55) 

Tbe physical masses fee the q. r{ and their LO mixing angle 0 from the fit are found to be 

ro n = (496.4 ± 1.3) MeV, m* “ (969.8 ± 5.8) MeV. 9 m -18.9° ± 0^=. (56) 

Tbe resulting plots can be seen in Fig. 1. We verify that if the physical masses are excluded in the fit. Af 0 - 
813 ± 11 MeV results. If wr only include tbe physical masses and exclude tbe lattice simulation data in the fit. 
Afo — 859 ±11 MeV b obtained. These determinations of Af 0 lie within the broad range summarized in Ref. (43) and 
are quite close with the commonly used values of Afo — 850 MeV |43|. Taking into account tbe large uncertainties 
of the lattice simulation data, specially for m^. and the concise formalbm of tbe LO mixing, it is impressive that 
the lattice simulation data can already be qualitatively described with the LO analysb, as shown in Fig. 1. This also 
indicates that tbe higher order mixing effects can only give moderate corrections to the masses of rj and q'. 
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NBWtkka, in order to describe the Uttice data more accurately, specially the ff masses. the chiral corrections 
beyond the leading order ore Deeded. For the physical masses, it has also been shown that the LO description fails 
to explain the mat* ratio of i) and if accurately enough |41]. Therefore it is essential to generalize the discussions to 
NLO and NNLO in order to achieve a precise description both fee lattice simulations and physical data. 
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FIG 1. TV mams of r? az*J if' ham the LO hi. TV left most two petal* correspond to the physical i u au« The remaining 
lattice simuktlcn data arc taken from Kali. |12. 13] (iTTOC), Jll] (UKQCD), [loj (RBC/UKQCD). |9| (USC). abac w* caly 
take into account the sliuilaUaa pouts with < 500 McV The shale area surrounding each curse Hands for the statistical 
uncertainty trees the fit. 


B. Next-to-leading order analyses 

At next-to leading order, in addition to the parameter \4q at leading order, there are five additional free 
the decay constant F at chiral and large Sc limits, and the four NLO LECa L*. X*. Ai and Aj in Eq. 
order, os well os at next-to-next-toleading order, ooe can rewrite the chiral expansico of the chscrvoblcs in vanoias 
equivalent ways up to the perturbative order in A under consideration. In the following disc us roan, we will perform two 
type* c£ fits: one usi n g F in the theoretical NLO and NNLO expressions and the other employing F n . os discussed in 
Sect. IID. Since the differences of the theoretical expressions used m the two types of fits are beyond the considered 
precision, the tnnanccs of the outputs from the two fits cun be considered as systematic errors from the theoretical 
models by neglecting higher order contributions. In the following, we will explicitly present the fit results by using 
F in the theoretical expressions, which is the moat straightforward option, as discussed in Sect. IID. The outputs 
of the fits with the theoretical formulas expressed in terms of F n will be uoed to estimate the systematic errors the 
difference bctwwn the central values of the two types of fits will be i»«l to estimate the truncation uncertainty due 
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to wcrking Just up to a given order in the A expansion, providing the seccod error for each quantity in the following 
tables. 

In Refs. [34-36). it is argued that at each chiral order, the leading Nq effects are dominant, or in other words that 
the Ai and Aa terms ore assumed to be much !e« irrelevant than the Z.& and L% terms in the NLO 6 expansion. This 
assumption has been mccc or lens confirmed when focusing on the masses of rj. rf and the LO mixing angle at the 
physical points (34-36j. In Ref [37], the local higher order LECs were estimated by the tree-level resonance exchanges 
and it was found that with those LECs A 2 5c< * ln3 to be more important than A 2 when focusing on the physical moni es 
for rj and rf. It is interesting to check how these assumptions work when including the lattice simulations and the 
phenomenological results of the two-mixing-angle parameters, which arc not considered in Refs. [34-37). Different 
sets of fits to the lattice data and phenomenological inputs freed the two-mixing-angle scheme ore performed either 
by fixing Ah-j.j to zero or releasing their values, in order to reexamine the assumptions. Interestingly we do not find 
qualitative changes between the fits with fixed Ai-jj = 0 and the ones with free values for these parameters. This 
teOs us that indeed the Aj and Aa terms do not significantly improve the fit results, even after taking into account the 
lattice simulaticcis. Nevertheless, we find that these two terms are quite important to reproduce the phenomenological 
mixing angles and 0* in the fits where A/o is fixed at its LO \ulue. If A/© is released in the fits wr find that including 
Aj and A 2 improves the descriptions of from Lattice simulations. Therefore, we will not further discuss fits with 
Aj and A 3 set to zero in the following. Instead, we focus on the results given in Tahle I with all the four NLO LECs 
in the fits, namely L t . I*. A 2 and Aj “ E<3- W- 

Fee the parameter A/o, wo take two strategies to estimate its value in NIX) analysis. In one of them we fix 
A/ 0 - 835.7 MeV from its LO determination (NLOFit-A) and in the other case wo free its value for the NLO fit 
(NLOFit-B). These two NLO fits are given in Table I. The first error bar for each fitted parameter corresponds to the 
statistical ooe from the fits and the second errcc bar is estimated freen the variation of the fits between those using F 
and F r in the NLO (and later also NNLO) theoretical expressions. From the two fits shown in Table I. one con see 
that releasing A/o in the fits barely changes the fit quality with respect to the cases when its value is fixed, although 
there are slight variations in the determinations of A/o and Aj. 

Concerning the results of the LECs in Table I. the resulting values for F from the two fits are quite compatible and 
close to the physical pion decay* constant. For A 2 and A 2 . their values are poorly kzxrw*n in literature and it is helpful 
to compare our values with the following estimate for their ranges: we take the LO determination A/o — 835.7 MeV. 
and we then separately include the A 2 and A 2 terms in the ry- rf mixing and vary thdr values to obtain new results 
for m, and f7j„- with the physical m*. Since the Aj and A 3 terms are NLO l/Nc effects, it is reasonable to assume 
that thdr correct ions to mj or m*. should be at most around 30% of the LO results. In this way we can set up 
conservative and rough estimates for the ranges erf At and A 3 . winch are found to be 


|A|| < 0.4. |A 3 | < 0.7 . (57) 

The resulting magnitudes of Aj in our fits are tiny and consistent with zero, as shown in Table I. For the parameter 
A 2 . our determinations lie within the ranges estimated in Eq. (57). Its value, spedolly the ooe from NLOFit-A, is dose 
to the one used in Ref. (48). where the mixing was discussed at next-to-leading order. However the determinations fee 
A 2 in Table I become much more precise than those given in Refs. [37. 39). where the lattice simulations for m p and my 
ore not included, indicating the usefulness of incorporating the lattice data in the £/(3) *PT study. Our determinations 
erf L* and Is are in goed agreement with the leading Nc predictions freen resonance chiral theory [49(. the 5t/(3) 
cue loop results in Ret |1) and the one-loop resonance chiral theory determination for Ls (50). But the \ulues here 
arc clearly larger than thosse from the recent two-loop determinations (3). 21 ). the results from Kj scattering in the 
scalar channels |51|. and the ooe-loop resonance chiral theory estimates for La [23). The discrepancies of La and Lj. 
comparing with the recent twx>-loop determinations (20. 21). can be eliminated once the Otff) LECs arc taken into 
account, as we will show in the NNLO discussion. 

The values of the parameters m the two-mixing-angle scheme and the mass ratio of strange and up/down quarks 
resmhing from the fits are given in Table II. Similarly, the first errcc bar for each quantity* is the statistical error and 
the second one corresponds to the systematic error, which is obtained in the same way as the one in Table I. Notice 
that these inputs have already been satisfactorily reproduced in NLO Analyses. 

The other quantities in the fits arc presented in Figs. 2, 3. 4 and 5. together with the lattice simulation data and 
the experimental inputs. We find that the final outputs from NLOFit-A and NLOFit-B ore quite similar, so only 
the plots from NLOFit-B arc given explicitly. The shaded area surrounding each curve car responds to the statistical 
error band for each quantity. In Fig. 2. we show* the resulting figures from NLOFit-B for the masses of rj and rf. In 
Figs. 3, 4 and 5. we show the corresponding plots for mj^. F v .k and F*/F* as functions of mj. respectively. 
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TABLE l Paranatal from the NLO dt* The moaning of different nototxns to label different fit* arc explained in detail In 
the tat. In the row of Af*. the column* with «».?• deoxc the fit result* by fixing the value U M 0 from it* LO determlnatvn 
The first cncc bar for cxch parameter 1 * the statistic*! one give* by the fit* and the second one correspond* to the systematic 
error The way to estimate the systematic error * explained in detail in the text. 



TABLE II The output* from NLO fit*. Notice that and if. arc not the pbcr*amcr»afotfc*l Inputs in the fit*, since 

they are related to ft. F*.4fc and through Eq. (24). The phcficmec^c 0 coi values for the m ixi n g parameters arc take* free* 
Red. | 1 !>] and we tuple the error bands here m order to make a conservative estimate The Input of m«/M to taken from the 
FLAG weekmg group in Red. |46| and w* ass^n the lirX error bar as done In Red [ 20 ]. ft* the mce bars of each quantity, the 
first one corresponds to the statistic crroc and the second one is fee the system*’.* error, which ore explained m derail In the 
text. 


C. NLO fits focusing on the masses 

In this section, we present another kind of NLO fits by focusing on the masses of ry.f f,K and excluding the decay 
constants F r . F* and their ratio. This kind of discussion is well motivated, since it is known that the NNLO 
corrections in 6 counting, such as the LEX? I 4 . ore important to simultaneously describe F r and F* [20—22]. But this 
LEC is absent in NIX) study. We have also provided another independent confirmation cq this finding in Fig. 4, where 
one can see that the decay constants cf pion and koon are poorly reproduced at next-to lead i ng order in 6 expansion. 
When only focusing on the q. rf and K masses and the ratio m,/m at next-to leading order the parameter F can not 
be resolved, because it always appears in the feem Lt/F 2 or L%/F 2 . We will fix its value to F = 90 MeV. close to 
the values given in Table I. For the mixing parameters we consider the mixing angles of 0 o and #*. but exclude the 
constants F 0 and F # . This is because F 0 and F, are dependent on the parameter F and should be determined together 
with Fr and F*. For simplicity in later discussion, we call the fits performed in this section os the mass-focusing type 
throughout. 

As in the previous section, we present the fits with F in the denominators of the theoretical expressions (c.g. 
Eq. (27)) and use the fits with F r to estimate the systematic errors. For each case, we perform the fits cither by 
fixing A/u at its LO determination (NLOFit-C) or by freeing its value (NLOFit-D). The fitted parameters are given 
in Tkble III and the m./th ratio and mixing angles are given in Table IV. The resulting figures from NLOFH-C and 
NLOFit-D are quite similar and we explicitly show one set of them. e g. NLOFit-D in Figs. 2 and 3 for the and 
kaon masse*, respectively. 

A significant difference between the results in Table I and the mass-focusing fits in Table III is that much larger 
statistical error bars are obtained in the latter case, especially for the LECs Lr, and La. as they are constrained by 
fewer data. Likewise, there are large systematic errors for the values of L*. and La in Table III. indicating a larger 
truncation uncertainty* due to higher orders. We do not see a significant improvement when freeing the value of Mo 
in the fits. 
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FIG. 2: The c^tu *i erf ) and if ***** the NLO and NNLO flu. The kft era* tw: pant* correspond to the physical mo&«* 
The rcmidiilc* lattice simulation data we taken flam Rcto (12, Ul (HTUC), [11] (UKQOO), (10) (RBC-UKQCD), (HSC). 
where wc only take into tc count the point* with m« < KIJ Me V. The shaded wcm wound the block solid and red doibed 
lines stand for the statistical moc bands horn the NLOFtt-B and NNLOFlt-B fits, respectively. The meanix* of notation* for 
dlllueut Um» we explained In detail In the teat 
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FIG. 3: Kaon ootf from the NLO aid NNLO fit* The lattice slxnulatxn data Are token from RBC and UKQCD (16, 17] 
Only the unitwy point* siiuuUted with the physical itionge quark uiua tit included The shaded area* orcailid the black sedd 
and red daibcd lines at sad lor the statistical csroc hoods from the NLOFtt-B and NNLOFtt-B fits. respectlvnly. The meanin* 
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TABLE IV The outputs from the dim*- focusing NLO fits See Tnbk U fee the phcMne&>k^Ksl inputs. The first error 
ice each quantity corrcsncmds to the statistical cue and the second csrcc denotes the systematic uncertainty. See the text far 
details. 


D. Next-to-next-to-loading order analyses 


From the NLO discussions m the previous two sections. we observe that the phenomenological results and the 
lattice simulations ca q and rf states can be rearonahh* reproduced. This is an important improvement comparing 
with the LO study, since at this order we only haw the conventional coc-mixing• angle formalism. The tw»mixing* 
angle formulism only shows up beyaad LO. However. observing m*. F». Fk atxl their ratio in Figs. 3. 4 and 5, it is 
clear that the NLO analysis is still inadequate. We need to include higher order contributions beyond NLO in order 
to further improve the descriptions. Moreover, the chiral logarithms predicted by xFT at one loop start at NNLO 























PIG. 4 Ptan and loan decay cowtanLs free* the NLO and NNLO flu. TV left-most points far F m and fW cowpcmd to the 
physical experimental input* IV remains* lattice slxanlatian data are taken flow RBC and l/KQOD [16. 17], where wo have 
c«Jy mcliaicd the unitary p<4nt* simulated with the physical stiange quark mas*. The ahadtd area around each curve stands 
Ice the statistical error band from the flu. TV means* of notations (or different Una are explained m detail in the ted 


in the 6 cxpuisy.it! Due to their importance xn ether observables. we consider it is relevant to ihstruso the impact of 
these chiral logs. 

As in the NLO cane, wr perform two types of fits, using the NLO and NNLO theoretical expressions given in terms 
of F and F w far %nrious observables Wr explicitly present the fit results with F in the theoretical rxprwuionn and 
use the alternative fits exprexed in terms of F r to estimate the systematic errors, due to war long up to NNLO in 
<5 and neglecting higher orders. According to the Lagracginn in Eq. (5). eleven additional unknown LECs appear 
at NNLO and there will be seventeen parameters in total for the NNLO study. At the present precision of the 
lattice simulations and phenomenological inputs, it is impossible to obtain sensible and stable fits if we free all nf the 
seventeen parameters Therefore, we need to take other independent determinations for scene of the LECs in order 
to proceed the NNLO study. 

\Ve mention thnt the stntc-ofart determinations of the 0(p 4 ) LECs in St/(3( *PT suffer uncertainties from the 
many poorly known Oijfi) LECs [20. 21). Because of the Large cumber of harely known C?(p c ) LECs. it is rather 
difficult to get conclusive r«iuits in the present two-loop 5I/(3) *PT studies [TO, 2I|. In the present work, there 
are five Oijfi) LECs. Le. Cu.Cu.C jt,Cu, Cm- in Eq. (5) and we cannot make precise determinations of these C» 
parameters here. Maybe when taking into account the scattering data, one nan make more stringent constraints on 
the C, LECs in t/(3) xPT. But this is beyond the scope of current work. Instead we take the C, values from the 
Dyson Schwinger like approach given m Ref. [52]. where ail c£ the 0{j> c ) C, at leading S c are predicted. In order 
to show the dependences of the final results on the C t tables. we alro perform other fits by using their updated 
determinations [53]. Like in Ref. [20]. wr multiply the 0(p c ) C t from Refs. |52. 53) by a global factor a and consider 
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cnee band from the fits. The meaning of notations for different line* are cxplalnoi m detail in the ton 


a *0 n free parameter in the fits. In this way. we partially compensate the large uncertain tie* of the C t parameters. 

For the operators proportional to i^ 21 . L is and L*\ in E«f (5). they are not present in SL T (3) *PT and purely 
contribute to the r^rf mi png, being irrelevant to the pion and kaoo observables. Since the iff/ mixing parameters 
have already been satisfactorily described in the NLO fits, we do not further include t»J J \ Li§ and ot NNLO study. 
2 Their indusico in the present analysis tend to make the fit unstable. Clearly studying more related observables 
it would be pce&sible to extract these parameters but this is out of the reach of the present analysis. A global fit is too 
unconstrained, being unstable and producing values o i the latter couplings compatible with xcro within uncertainties. 
Then we are left with three 0(N%,p 4 ) operators. L 4. I« and L? . which have corresponding parts in 5L r (S) yPT. Since 
1/(3) and SL r (3) X’FT ccertain different dynamical degrees of freedom, the corresponding LECs from the two theonts 
can be different. A typical example is the Lj parameter in Slf(3) \PT, which is demonstrated to be do mina ted by 
the singlet »» state |1|. Since in f/(3) xPT the singlet r \j has been expbdtly introduced, the value oc Ly in this theory 
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con be totally different frccn I^ 13 ’ in 5t/(3) cose. While for other 0 {p A ) LEcs. such os Ltw4.tA*, the diffexencre 
between 6/(3) and 56/(3) \FT ore not expected to be os large os the Lr case. since they do not receive the tree-level 
contributions from the o© state. 

Another subtlety to take into account is that m, and appear in the chiral lcope and, at the ^ne time, the 
final expressions of m„ and m,» depend on the three loops os well. In ccdcr to avoid making the complicated iterative 
procedure to obtain the ryif mixing parameters, we use the LO formulas fee m n and m,. in the chiral loops. The 
differences caused by this simple treatment and the strict iterative procedure arc beyond the NNLO precision in 6 
expansion, since the chiral loops themselves are already NNLO. Our simple sdution is also justified by the fact that 
the IX) description of m v and is in qualitative agreement with the lattice simulation data, os shown in Sect. Ill A. 
Since the qualitative agreement between the LO formulas and the lattice Simula tico data requires the value of M 0 
to be around 835.7 MeV. as given in Eq. (55). we fix A/j - 835.7 MeV in the following discussions. This also helps 
to stabilise the NNLO fits, with its many free parameters. Other useful criteria to discriminate reasonable fits arc 
the a priori ranges estimated in Eq. (57), since the fits with large magnitudes of Aj and Aj imply imphysically Urge 
corrections to the iyi{ mixing parameters and the breakdown of the 6 expansion. In the following we only present 
the fit results that are consistent with Eq. (57). 

With all of the above setups, the values of parameters from the NNLO fits are summarized in Table V. The fits 
labeled by NNLOFit-A and NNLOFit-B correspond to using different values of the Oip *) LECs. For NNLOFit-A. 
the C values are taken from Ret [52]: 

Ci2 - -0.34. Cu m -0.83. CiT - 0.01, C J9 = -0.48, Cn - -0.63, (58) 

which are given in units efi 10” 3 GcV“ 2 . Rir NNLOFit-B. we take their updated 0{p 6 ) C t valure from Ref. |53|: 

C U = -0.S4. C u m -0.87 1 C iy =0.17, C 19 =-0.27, C„ = -0.46, (59) 

in the same units as before. 

It is clear that the parameters resulting from fits with different C, inputs slightly differ from coe another. We 
remind that the first error bar for each parameter in Table V corresponds to the statistical ooe directly from the 
fits and the second error bar stands for the systematic ooe. which is estimated, as usual, from the variation of the 
parameter from the NNLO fits with the theoretical expressions in terms of F and those exprereed as functions of F r . 

At NNLO. one has the contributions from the chiral loops and the Oip*) LECs. which moke our determinations in 
Table V closer to the recent two-loop results of the 56/(3) *PT LECs, comparing with the NLO determinations in 
Table I. Some typical trends cf the values of parameters from the NLO study in Table I to the NNLO one in Table V 
are summarized now. The axial-vector decay constant F at leading Nc and chiral limit is reduced at NNLO. which 
is mainly due to the inclusion of L 4 . Our conclusion is based an the fact that strong correlations between F and L 4 
always appear, which has been confirmed in previctis study (22. 23]. Fee L\ and L%, we find that their values arc 
obviously reduced compared to the NIX) determination and become closer to the two-loop results in Ref. (21]. As 
mentioned in the former reference, the discussions in the two loop 56/(3) \PT are sensitive to the value of the \/Sc 
supprereed LEG L A . The present study provides on independent determination for this parameter and for the \/S c 
suppressed LEG Lt os well. We mention that our determinations cf L 4 have opposite signs with respect to that in 
Ref. (21|. which may be the source of the smaller F obtained in that reference. Notice that the present valure of 
L 4 , Lt,, Lt, L% are rather compatible with the combinations of 2L§ - L% and 2 L t - La given in Ref. |54|. Fit solution 
with larger As and Aj than those in Eq. (57) (out of the a priori range (57)) are discarded: they are not considered 
as reasonable physical solutions and will not be discussed any further. According to the values of a in the two fits, it 
seems that our study somewhat prefers smaller magnitudes of the 0(p l ) LECs than these from the Dyson-Schwinger 
approach given in Refir (52, 53) and also prefers a global change of sign with respect to Eqs. (58) and (59). We have 
investigated the impact of fitting a but releasing ooe of the 0(p c ) LECs as an independent parameter (c.g.. Cu), but 
no definitive conclusion could be extracted. These puzzles cannot be resolved here and it is definitely interesting and 
necessary to further investigate the valure of the 0(p 6 ) LECs in the future. 

The various plots from the NNLO fits arc shown in Fig. 2 foe m ? and m n >, Rg. 3 for m K . Fig. 4 for F r and F x , 
and Fig. 5 for the ratio Fx/F*. together with the NLO results and the lattice simulation data and experimental 
inputs. The shaded area surrounding each curve represents the statistical error band. The figures frccn NNLOFit-B 
are compatible with these from NNLOFit-A within the uncertainties, so we only show the results for the former in 
Figs. 2, 3. 4 and 5. 

In addition, to demonstrate the effects by using the LO masses in the higher order corrections, instead of the 
renormalized oore, we explicitly show the results for m, and m T expressed in terms of the IX) masses TTTr.TTTx and 
1 /F 2 in Fig. 2. with the lines labeled as NNIX)Fit-BTTr T *. The values of the LECs when plotting these lines arc 
exactly the same as those from the NNLOFit-B column in Table V. In this way, one con directly see the differences 
due to the N 3 LO truncation uncertainty caused by using the renormalized marere and the LO ones at the NNLO 
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level. According to Fig. 2, we conclude that the differences for m, and m v caused by using different types of maai w 
in the higher order corrections are rather within the statistical uncertainties freen the fits and therefore the differences 
should be perfectly compatible within the total uncertainties after taking into account the systematic ones in Table V. 
We verify that similar conclusions are obtained for other cases. In order not to overload the plots in other figures, we 
shall not explicitly show the results given in terms of 17T r and TTT*. 

From Figs. 2, 3. 4 and 5. we observe, when compared with the curves of the NLO study, slight improvements in the 
reproduction of the masses for ri.rf and significant ones for m K . F w , F K and the ratios of F K /F W . Moreover the \ 2 
for the NNLO fits are greatly reduced compared with \ 2 fee the NLO ones, indicating that the NNLO corrections are 
important at the present level of precision and essential to sinaaltaneously describe the lattice s i mul ation data and 
experimental inputs of the light pseudosealar mesons n. K, q and rf. 
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TABLE V: Parameters from the NNLO fits. In all * these fits. Af 0 to fixed at 835.7 MeV from ts LO dctermlnattco The 
meaning of different notatxns to label different fit* are explained m detail In the text. The first error bar Ice <och parameter 
cones pends to the statistical one and the second cnor denotes the systematic uncertainty See the text for details 
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TABLE VI The outputs from NNIX> fit*. See Table 11 for the cxplanaticn of the phenomenological inputs. The first error fox 
each quantity cones pends to the statistical one and the second cnor denotes the systematic one See the text for details 


IV. CONCLUSIONS 


In this article we have performed a thorough study an the r^rf mixing, and axial-vector decay constants for the pion 
<n, up to next-to-next*to-leading order in 6 expansion within C/(3) chiral perturbation theory. We have carried 
tailed scrutiny and discussions of our results, which have been carefully compared to other works in literature 
rf-rf mixing. A general mixing formalism, including the higher-derivative terms and kinematic mixing coses, 
n addressed in detail. The conneetions between the mixing parameters from the popular two mixing-angle 
and the low energy constants from chiral perturbation theory have been established, both for the singlet-octet 
d the quark-flavor basis. 

xjcsidered quantities, including the mass es of q. q r and K . the quark m v a ratio of m ( /m. the parameters in 
the two-mixing-angle scheme and the it. K decoy constants have been confronted with recent lattice simulations and 
phenomenological inputs. We find that the next* to leading-order fits yield satisfactory descriptions for the ma*&cs 
erf the three psetdeoealar mesons as functions of mj and the four mixing parameters (Fo,Ft,$o, 0»|- producing in 
addition reasonable values of low energy constants. Nonetheless, when the ir and K decay constants are included 
together with the masses and mixing parameters in the fits, the next*to-leading-order analyses are inadequate and it 






















is necessary to step into the next-to-next-to-leading-order study. Using tbe 0(p ‘*) LECs determinations from a Dyson- 
Schwinger-like approach [52, S3) multiplied by a global factor, we are able to achieve a reasonable description for all 
of the physical quantities considered above and tbe resulting values for the leading 0{p 4 ) low energy constants L, 
and £« turn to be compatible with tbe very recent two-loop determinations in Ref. [21). Therefore we cooclude that 
the large Nq 1/(3) chiral perturbation theory oSers a concise theoretical framework that is able to simultaneously 
reproduce accurately the general ty-r/ mixing and to provide sophisticated enough expressions to describe tbe chiral 
extrapolations of the r and K decay constants and m asses. 

Our results are also useful for future phenomenological studies of different processes involving q and rf. Combining 
Eq. (21) or Eq. (23) with Table VI. one can directly find the relations between tbe physical states i),q' and the 
octet-singlet bases rj*. »>3 or tbe quark-flavor bases q,. These relations are consistent with the requirements from 
the recent lattice simulations and phenomenology. 

Finally, it is worthy to remark that some of the parameters in our best analysis (NNLOFit-B) in Table V have been 
determined with relatively small errors. Rjr instance, the NLO parameters Ajj, which are fitted up to 0{N£ 2 ) in 
the NNLO analysis, become 


As = -0.04 ± 0.ft> ± 0.13 f A 2 ■ 0.14 ± 0.10 ± 0.40. (60) 

Tbe NNLO fit also determines some 1/(3) NNLO couplings with relatively high precision. NNLOFit-B yields 

10 * x L A = -0.09 ± 0.06 ± 0.23. 10 s x U = 0.03 ± 0.03 ± 0.02 f 10 3 x L 7 = 0.36 ± 0.06 ± 0.12. (61) 

Even though tbe error estimates in the present article must be considered with some caution, as some lattice systematic 
uncertainties escape our control, this hints the potentiality of this 1/(3) *PT framework. We hope these results may 
encourage future lattice analyses along this line. 
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Appendix A: Higher order corrections to the T? and Tf bilinear terms 


In the following we provide the explicit expressions of the £g's in Eq. (11). When expressing the results in terms of 
F, they take the form 
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~ 2 ^^I F 1 { ^< 8m * - 5mJ) ♦ 4(&">* + mi)«, * 3v^(-4nii; * m 2 )# 

+4f,(-Sm}f ♦ 2mi|»3 + 4V5(-mjf + m;)j*j -VmJ) 

«*&&•& - ml) + 3v^J(4mi - m 2 ,),‘, + <»(*"* + mj).’ + v^-8m^ + Sm^.J 

+44^-5-.*,+a-4)]-Mmj) 

♦6{V^J(2mJ f - mj) * c${im a K + mi)* + v7(-2m^ + mil^A^m},) 

♦9mi(-♦ c,, t ♦ v^»J)i4o(mi)J 
fl, /o»> .,(» 32Li(2m« c - - mj){%fej - c», - y'l'j) 

j2V^-imi(mi - mi) + - im‘ K m 2 . * 3mJ)i, + 2V£nJ f (-mJ f + mi)*2j 

2v / 7d(2mi - mlmj - mi) + c#(-4mi + 20mimJ - 7m*)s* + 


L6I 


lbL 7 
3^ 

2 v'I(- 2 m} ( + mj^m' + m*)»J 

512(1, - 2L,)Um 2 K (2m* K -m^m 2 - m*)(v/fcj - c,», - V2.J) 

37®- 

16 <I, - 2 U)A a [VfcJ(mJ c - mj) + c( 2 m' K + m‘) M + + mi),J| 

-37^- 

32{3C„ + 2C„)(4m‘ ( - 6m^mJ + 3 m^m* - mSHv^ - c,», - v7.J) 

3 ?® 

When expressing the above results in terms of F r from Eq. (27), the terms with L\ and L\L% can be different from 
the expressions in terms of F and the other parts remvin the same, apart from the obvious replacement of F by F r . 
Therefore, for the expressions of 6* expressed in F r . we only give the parts that arc different from those in terms of F 

= •^^■W(2m} i + 2m‘ K m 2 - m*) .2v&*«r("*t + m^m} - 2m*) + ^(mj, ♦ m^mj + m*)|, (A10) 


A9> 


C W - + mJcmi + m‘) - * m 2 K ml - 2m*.) + «J(2mt + 2m‘ K m 2 , - m*.)|. (A11) 

4r.).ti HBtjlmj, + mjml - 2mj)(^ - c,5, - vfoj) 

<*, =- m (Ai2> 


,<K>X.l. _ 123£,t, 

*< -3FT 


d(16m5f - 16m}fmJ +9mj) + 16v2c#m}f(mjf - mi)s* 
4(8mt - Smlm* + 9m£)s2 


[AW) 


•<F.).UL% 128L S I* 

m i 3Fi 


d(8mjf - Smjffni + 9m*) - - mJ)j* 

+(16mt - 16mtmi ♦ 9m?)»il . 


[A14> 



Ar.iX.im _ 10 2*LiL,m 2 K (m^ - )(J2c* -c r >,~ JU\) 

3f?- 


iA15> 


In ct<1« to obtain thr full 
precision required. Taking <5* for 


for the Vs given in terms of F r one has to make use of Eq. (27) up to the 
pie. its final expcemico in terms of F n » 


16 tt(m^ - miXvToj - »», - yfoj) 

flP - «»5#Aj 

C5»Au(mi) - mj!) + c(2mi + m})* + sffl-nA * mi)»il 


l-ir'/'i 


r 


- c.». - V%2) 128tj(m* ♦ m'mj - 2m*)(v^i - c,., - V%J) 


3f? 

64(C U + Cir)m^(m^ - m*)(v/^3 - c**, - V$»J) 

3f? 


If? 


(A16) 


which differs from Eq. (A6) in the L\ term. Ifor 6 \ . 62 and Jj. their expressions are the same regardless of whether F 
cc F r is chosen up to next-to-aext-to-leading order. 

Far completeness, we also give the results in terms of the LO masses 777 r and W K and i/F 2 . Only the terms with 
L x L y being L % and L } the NLO LECs in Eq. (4). will be different, comparing with the expressions in terms c£ m n 
and m/f and the other parts re ma i n the same, apart from the obvious replacement of the renormalized masses by the 
LO ones. Therefore, we only give the parts that are different from those in terms of m T . m*c and l/F 2 and it turns 
out that in this case all of the LtL } terms for vanish. 
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